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We consider the excitation spectrum and stability of quantized vortices in a weakly-interacting Bose gas trapped 
in a toroidal container, and discuss the driven rotation of such a condensate. 
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Quantization of circulation and the related phe- 
nomenon of persistent currents are two of the most strik- 
ing properties of superfluids. Q The recent creation of 
atomic Bose-Einstein condensates has renewed in- 
terest in the superfluidity of confined Bose systems. 
Here we discuss two related issues that pertain to the 
demonstration of these phenomena in trapped gases - 
vortex stability and the driven rotation of condensates. 

Harmonically trapped Bose gases (and, more gener- 
ally, gases confined to monotonically increasing poten- 
tials) cannot support stable vortices the absence of an 
externally imposed rotation. Rj In such a potential the 
vortex core is not pinned, and can slip out of the sys- 
tem, leaving behind a non-rotating condensate. |l0| The 
simplest trap that could pin a vortex core is a toroidal 
container, which confines the gas to a non-simply con- 
nected volume. For theoretical ease, we consider a dilute 
Bose gas in a narrow torus whose cross-sectional area A 
is small compared with the square of its radius R. This 
limit allows us to neglect transverse motion and to focus 
our attention on the more interesting azimuthal direc- 
tion. We regard the system as a ring Jll| of radius R, 
and label single-particle states by their azimuthal angular 
momentum. We use the Bogoliubov approximation p^j 
to explicitly calculate the properties of the dilute Bose 
gas in this geometry. 

Consider a collection of N bosons of mass M in- 
teracting via a (5-function pseudopotential of strength 
U = A-Kh 2 a/M, where a is the s-wave scattering length. 
The Hamiltonian is 
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where e q = hajq 2 is the energy of the single-particle 
state of angular momentum qh and V = 2nRA is the 
volume of the annulus. The characteristic frequency 
u> = h/(2MR 2 ) corresponds to the energy difference be- 
tween the ground and first excited state of a single par- 
ticle on the ring. 

First, we discuss the stability of a Bose gas prepared 
in a quantized vortex of circulation m(h/M), that is, a 
Bose condensate whose condensate wavefunction varies 
as e m ^. In the presence of a weak, static, asymmetric 
perturbation, angular momentum is not conserved. In 
a rotating normal fluid, scattering of particles against 



such an asymmetric perturbation will ultimately bring 
the fluid to rest. A characteristic feature of superflu- 
ids, however, is the absence of accessible final states into 

We show that unlike the 
, the pinning provided by 



which the vortex can decay. Q 
case of a monotonic potential 
the hole in the torus stabilizes the vortex if the scattering 
length is sufficiently large and positive. 

Second, we examine the production of a vortex start- 
ing from an initially non-rotating Bose condensate. We 
discuss the response of an m = condensate to a weak 
time-dependent potential that drives the system to rotate 
at angular velocity O. At small £1, a weak, adiabatic per- 
turbation imparts angular momentum to the condensate, 
but the corresponding moment of inertia / = (L z )/Q is 
proportional to the square of the perturbation strength 
and can be small compared with /classical = NMR 2 . This 
phenomenon of non-classical moment of inertia Q is es- 
sentially the "Meissner effect" for superfluids, which "ex- 
pel" rotation in the same way that superconductors expel 
magnetic fields. We determine the critical angular veloc- 
ity in the toroidal geometry, which is simply related to 
the Landau critical velocity, and discuss the creation of 
vortices when this critical velocity is exceeded. 

Excitation spectrum of a vortex. We begin by de- 
termining the excitation spectrum of a quantized vortex 
in a dilute Bose gas, in the lab frame. Although this 
spectrum can be derived by transforming to and from 
the co-rotating frame of reference, JHJ it is instructive 
to calculate it directly in the lab frame using the Bogoli- 
ubov approximation. []l2| , [l5| Following the familiar pro- 
tocol, we replace the creation and annihilation operators 
for the vortex condensate by a c-number, and retain only 
quadratic terms in the remaining creation and annihila- 
tion operators in eq. (0). On a ring the vortex condensate 
of quantum number m is the unique state e m8 . |l6| The 
Bogoliubov Hamiltonian becomes 
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where Ec{m) = e m + Un/2 is the Hartree-Fock en- 
ergy per particle, n = N/V is the number density, and 
we have substituted the appropriate chemical potential 

Mm "T" UTl. 
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The Hamiltonian (0) is easily solved by a canonical 
transformation to the quasiparticle operators 
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where the amplitudes u q and v q must satisfy the normal- 
ization condition \u q \ 2 ~ |w 9 | 2 = +1 to ensure that Bose 
commutation relations [/3 g ,/3|] = 1 are preserved. Note 
that excitations are labeled by the angular momentum 
relative to the condensate, and are not states of definite 
angular momentum in the lab frame, Jlj|[| due to anoma- 
lous scattering against the rotating condensate. There is 
a single quasiparticle excitation for each non-zero (i.e., 
positive and negative) integer q. 

such that /3| is an 



We seek the amplitudes u 



eigenoperator of 7iBogo with excitation energy hu q , i.e., 
[H-Bogo , fi q ] = hk>q0 q - This leads to the (non-Hermitian) 
eigenvalue equation 
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which has the two eigenvalues 



^ ± y/e(m, q) 2 + 2e(m, q)UN, (5) 



by the superflow. If the q = ±1 modes are excited co- 
herently, the resulting density fluctuation will precess at 
angular velocity — w_i]/2 = 2uro, which is reminis- 
cent of Vinen's original detection of quantized vortices. 

Persistent currents. We may use the spectrum (^) 
to assess the stability of a vortex of quantum number m 
in the presence of a static, weakly asymmetric perturba- 
tion. We assume that the gas has been prepared in a 
vortex state at low temperatures. A necessary condition 
for persistent currents is that the excitation spectrum of 
the current-carrying state must be positive for all relative 
angular momenta q: otherwise, even a weak asymmetric 
perturbation could excite a macroscopic number of such 
excitations, and the condensate would be destroyed. In 
the non-interacting limit Un — ► 0, vortices are always 
unstable in any geometry, since the total energy of the 
system can always be reduced by transferring particles 
from the rotating condensate to the (non-rotating) sin- 
gle particle ground state. Thus vortex stability depends 
critically on interparticle interactions. 

A vortex of angular momentum m > will be stable if 
u q is positive for all q, which ensures that excitations will 
not proliferate. It is easy to see that if u> q= -\ is positive, 
then all other tu q will also be positive. A little algebra 
then establishes the stability condition 
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where e(m, q) = \[e m +q + Cm-q — 2e m ]. For particles on 
a ring (but not, for example, in a harmonic potential), 
e q = hu>q 2 , and e(m, q) = e q . 

Only the positive sign in eq. (j^) permits the normaliza- 
tion \u q \ 2 — \ v q \ 2 = +1. The eigenvector corresponding to 
the negative sign has \v q \ > \u q \, and cannot be normal- 
ized in this manner; it therefore cannot represent a quasi- 
particle excitation. [fl7| This unnormalizable solution in- 
stead represents the quasi/ioZe operator, /3_ 9 , which is the 
Hermitian conjugate of the (normalizable) operator f3_ q 
that creates a quasiparticle of relative angular momen- 
tum — qh. From eq. (jij), u~ q = u q , and v- q = v q ; for a 
ring, these coefficients are independent of m. 

We conclude that in the lab frame the energy of a 
quasiparticle of relative angular momentum qh above a 
vortex of quantum number m is 
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^2Une q + {e q ) 2 . 
In the strongly interacting limit Un/huj 3> 1, 
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huj q m ^ m 2hLumq + V2Unhuj\q\ + 0{q 2 ). (7) 

The frequencies uj+ q and are different in a rotating 
condensate, reflecting the broken time-reversal and par- 
ity symmetries of the vortex. This implies a "Sagnac 
effect" [jl8| - if we create clockwise and counterclockwise 
propagating density fluctuations, they will move at dif- 
ferent velocities in the lab frame as they are carried along 



Evidently the dimensionless ratio of the mean-field inter- 
action to the level spacing must be sufficiently large for 
vortices to be stable. This requirement reflects the fact 
that along the decay path of the vortex there must be an 
intermediate state with a node along the ring at which 
the winding number of the phase "slips" by 2tt, and the 
vortex "core" passes out of the torus. M The energetic 
cost of introducing this density fluctuation depends on 
the compressibility of the gas; for sufficiently strong re- 
pulsion there is a barrier to the phase slip, and the vortex 
becomes stable. 

Negative scattering length. It is interesting to con- 
sider the case of negative scattering length, as is appropri- 
ate for 7 Li, where Bose condensation of a limited number 
of atoms in confined geometries has been demonstrated. 
H The in = condensate is stable with respect to quasi- 
particle excitations as long as the argument of the square 
root in (^|) is positive, i.e., when interactions are suffi- 
ciently weak compared with the level spacing: 
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This is simply eq. (g), with m — and U < 0. When 
eq. (||) is violated, uj q becomes complex, signaling an ex- 
ponentially unstable mode. A similar requirement holds 
for the harmonically trapped gas |l9| ]. 

Similarly, we find that gases with negative scattering 
length in a toroidal geometry can never support persis- 
tent currents, since vortices are then always unstable: 



2 



u)g—-% is always either negative or complex when U < 0. 
Again, this result follows from eq. d§|), with m 7^ 0. 

Making vortices. A vortex state can be created by 
cooling a Bose gas through its condensation temperature 
while steadily "stirring" it at sufficiently large angular ve- 
locity Q, which can be accomplished with a rotating per- 
turbation V(9 — fit). In the limit of a weak perturbation, 
single-particle states can still be adequately labeled by 
angular momentum. In a steadily driven state the macro- 
scopically occupied single-particle wavefunction will be 
that which minimizes the single-particle energy in the 
co-rotating frame, e' m — e m — Mlm. For example, when 
lu < f2 < 3cj the m = 1 vortex will be formed at the 
Bose condensation temperature. When the rotating per- 
turbation is turned off, a vortex in a toroidal trap will be 
stable if eq. (||) is satisfied. 

Landau criterion. Can a vortex be generated by 
stirring an initially non-rotating condensate? To simplify 
the discussion we assume that the perturbation is weak, 
and consider initial temperatures well below T c so that 
we may continue to use the Bogoliubov approximation. 
The time-dependent rotating perturbation becomes 
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NV q (u q + v g )e- in ^(pl+f3. g ), (10) 



where V q is the Fourier transform of the perturbation 
V{9), and we take Q > 0. 

There are two simple limits to consider. First, we may 
consider turning on the rotating perturbation adiabat- 
ically, i.e., V(6,t) = [1 - e~ t / T ]V(9 - tit), where the 
"turning on" time r is long compared with \u> — fi| — . 
Inverting the Bogoliubov transformation (|^) , we see that 
the condensate becomes 

mt) = (m) = -^{(b^+Ys^M (n) 

V2irN q ^ Q 
= {VN + Y,e lqe [u q (P q )+v q (Pl)}}, (12) 



where by first order time-dependent perturbation theory 



= VNV q (u q + v q ) c _ iQqt 
h(uj q — Clq) 



(13) 



when the perturbation is fully developed (i.e., t^$> r). 

The condensate responds at the driving frequency f2 
by forming density minima and maxima that track the 
peaks and valleys of the rotating perturbation. (This 
is easily seen by considering "J (6* — fit), i.e., in the co- 
rotating frame.) If V(6) = 2V q cos(qd), then the steady 
state ty(6,t 3> t) has non-zero angular momentum 
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but 4" only describes a vortex when its phase winds 
around by 2tt, which requires that the coefficient of e 10 



be large compared to that of 1 and e~ lB . For a weak 
perturbation this requires the near resonance condition 
il ~ Lo q /q. The energy transferred to the gas per unit 
angular momentum, [(H) — Eq]/ (L z ), is (lu 2 + £l 2 )/2qf2, 
which tends to u> q /q near resonance. 

If the same perturbation is turned on suddenly, it in- 
duces transitions between states that differ by energy 
TiVLq and by angular momentum hq. This resonance con- 
dition requires Q — to q / ' q for some q. For the Bogoliubov 
spectrum (@) with m = 0, the minimum value of u> q /q 
occurs for q = +1, and the critical angular velocity is 



critical 
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Un 



(15) 



In the limit of large Un/hu, ft crit icai -> [Un/ '(MR 2 )] 1 / 2 . 

Eq. (|l^) is nothing but the famous Landau criterion 
Jl3|] applied to the dilute gas in a rotating torus. For a 
dilute gas, there is no roton minimum, and the Landau 
critical velocity is the speed of sound c = [Un/Ad] 1 / 2 . 
The velocity v por t of the rotating perturbation relative 
to the (non-rotating) condensate is £IR, so according to 
Landau's argument the driven Bose system will be stable 
with respect to the proliferation of elementary excitations 
as long as QR does not exceed c, which recovers (|l^) in 
the large Un/fuv limit. 

The driven gas. When the gas is driven at or 
near resonance, its angular momentum and energy both 
grow beyond the linear response regime. To understand 
the fate of the driven gas, it is useful to first consider 
the (many-body) ground state of a Bose gas subject 
to the constraint that its total angular momentum is 
M z h. The overall ground state is always at M z = 0, 
since the true ground state must be nodeless; for small 
M z where quasiparticle interactions may be neglected, 
E G (M Z ) = E G (0) + hu q=1 \M z \. 

For sufficiently large Un/tko, the ground state energy 
Eq(M z ) will have more than one local minimum. A full 
calculation of Eg(M z ) is beyond the scope of this paper, 
but we may understand its general features by consider- 
ing the family of condensates of the form 



* x {0) - — WT=x- 

Ztt 
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where < x < 1. The total angular momentum of ( |l6| ) 
is (M z ) = Nx, and its total (Gross-Pitaevskii) energy is 
E G (M Z ) = htuM z +\Un[N+2M z -2M 2 /N], which devel- 
ops a second local minimum at M z — N for Un/hu > 1. 
The barrier between the minima at M z = and N arises 
from the presence of a node in ^ when the winding num- 
ber of t he p hase changes by 2tt, which occurs at x = 1/2 
in eq. (|lq). The density fluctuation that accompanies 
the node at intermediate x costs energy due to the com- 
pressibility of the interacting gas; for small and large x, 
however, the magnitude of ^> is nearly uniform, and no 
additional energy costs arise. Since d 2 E G /0M 2 < for 
Q < M z < N, only the M z = or N states can be 
thcrmodynamically stable. 
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Armed with this understanding of Eq(M z ) we can de- 
velop a scenario for spinning a stationary condensate 
into a vortex. First, we imagine an adiabatic pertur- 
bation which increases the angular momentum of the 
gas while keeping it in the corresponding ground state. 
Then during a short time interval the energy increment 
SE and the angular momentum increment 5M Z must 
be related by SE — [dEc / dM z ]5M z , which requires 
that the angular frequency of the perturbation vary as 
Q(i) ee [dE G /dM z ]\ MzW . (The absolute values of SM Z 
and SE depend on the strength of the perturbation.) 
Since Eq(M z ) has a maximum, the angular velocity f2 
diminishes as the total angular momentum increases. For 
our simple model, Q(M Z ) =uj + Un~- 2Un(M z /N). 

Eventually, the maximum of Ec(M z ) is reached. If the 
angular momentum is increased further (by continuing to 
drive the gas at low frequency), the gas enters a regime 
in which BEq j dM z becomes negative, and the ground 
state energy decreases with further increase in M z . Since 
the gas is thermodynamically unstable, in the absence of 
any further drive it will relax towards the nearest free en- 
ergy minimum, which is now the vortex condensate with 
m = 1. The energy liberated in this relaxation process 
goes into generating excitations above the ground state, 
and heating this normal fluid. The angular momentum 
gained by the condensate appears at the expense of an 
opposite angular momentum change in the normal fluid, 
which ultimately decays due to scattering by small inho- 
mogeneities in the trap. 

A back-of-the-envelope estimate suggests that heating 
due to non-adiabatic driving of the condensate is not 
a problem. If we simply drive the gas at frequency ~ 
^critical long enough for it to acquire angular momentum 
Nh, then it will absorb an energy of order ~ AT7if2criticai- 
At this point we rapidly turn off the perturbation (or turn 
its frequency down to uj < Q < 3w as discussed above), 
and allow the gas to reach a new equilibrium. If the gas 
has not heated up past T c , then it should relax into the 
m = 1 vortex. The condition is that the total energy 
transferred during the excitation period, ~ Afi,f2 cr iticai, 
must be less than ~ Nk B T c + E G (M Z = N). Taking 
ksTc ~ h 2 n 2 / 3 /M as appropriate for a dilute gas of den- 
sity n, we find that a vortex can be created in this manner 
if n > [a 3 /R 6 ]. If R ^> a, there is no conflict with the 
dilute gas condition na 3 <C 1. 
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E GP m = J dr[**# *+ ^|*| 4 ]. (17) 

Consider the energy of a condensate \& + <j>, where 
SEgp/S^ = [ity and <f> is a small (complex) variation. 
To quadratic order in <f>, the variation in Eqp is 



dr 



[H - v]<f> + ^[* 2 (^*) 2 + (**)V + 4|*| 2 |0| 2 

(18) 



For a vortex on a ring, = e lm6 , and this quadratic form 
can be partially diagonalized by defining 

4> q {6) = a,e l(m+9)s + &*e i(ro -« )8 (19) 

where (a q ,b q ) is the eigenvector of the (Hermitian) ma- 
trix 



- e m + Un] Un 
Un + Un] 



(20) 



Note that ( |19[ ) and ( J2C| ) are related to but distinct from 
(^|) and Q). The vortex is then locally stable if the two 
eigenvalues 

5 = e q + Un ± [(2m%) 2 + (Unf] 1/2 (21) 

are both positive for all q. This condition is the same 
as (^). More generally, for arbitrary traps one finds that 
the positivity of the elementary excitation spectrum co- 
incides with the stability of the Gross-Pitaevskii energy. 
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